Journal  of  Power  Sources  222  (2013)  359-366 


ELSEVIER 


Contents  lists  available  at  SciVerse  ScienceDirect 

Journal  of  Power  Sources 

journal  homepage:  www.elsevier.com/locate/jpowsour 


P  r*! 

Sbb..«ltS 


Elastic  configurations  of  self-supported  oxide  membranes  for  fuel  cells 

K.  Kerman3  *'* 1 *,  T.  Tallinen3’cl,  S.  Ramanathan3,  L.  Mahadevan3'b 

a  School  of  Engineering  and  Applied  Sciences,  Harvard  University,  Cambridge,  MA,  USA 
b  Department  of  Physics,  Harvard  University,  Cambridge,  MA,  USA 
c  Department  of  Physics,  University  of  Jyvaskyla,  P.O.  Box  35,  FI-40014  fyvaskyla,  Finland 


HIGHLIGHTS 


►  Address  the  origin  for  far-from  threshold  buckling  patterns  in  self-supported  oxide  membranes. 

►  Provide  an  outline  of  shape  evolution  based  on  reduced  stress. 

►  Corroborate  results  of  elastic  model  with  experimental  fabrication  of  self-supported  thin  film  YSZ  membrane. 

►  Utilize  self-supported  membrane  to  show  a  proof-of-concept  thin  film  SOFC  operating  directly  with  natural  gas. 

►  Obtain  near-ideal  open  circuit  potential,  indicating  pristine  nature  of  membrane. 
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Ultra-thin  oxide  films  are  of  interest  in  energy  conversion  technologies  such  as  low  temperature  solid 
oxide  fuel  cells  and  permeation  membranes.  Understanding  their  thermo-mechanical  stability  is  an 
important  problem.  Edge  clamped,  self-supported  thin  film  membranes  show  hierarchical  wrinkles;  with 
the  largest  wavelengths  in  the  center,  while  smaller  ones  arise  near  the  clamped  boundary;  corre¬ 
spondingly  the  largest  strains,  with  tensile  stress  comparable  to  the  residual  stress,  are  in  the  vicinity  of 
the  clamped  boundary.  Our  results  can  be  understood  by  simple  scaling  arguments  and  are  valid  for 
membranes  in  the  post-buckling  regime  far  from  threshold.  We  confirm  the  validity  of  our  analysis  by 
quantitative  experimental  comparison  to  self-supported,  square  micro-machined  yttria-stabilized 
zirconia  membranes  of  edge  length  160  pm  fabricated  by  lithography.  The  modeling  and  experiments 
combined  provide  a  foundation  for  designing  failure  resistant  self-supported  membranes  of  interest  to 
energy  conversion.  We  show  this  by  utilizing  such  membranes  to  fabricate  thin  film  solid  oxide  fuel  cells 
and  demonstrate  power  generation  utilizing  natural  gas  as  fuel  at  ~400  °C. 

©  2012  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

Ultra-thin  membranes  are  of  broad  interest  in  electrochemical 
energy  devices  that  are  generally  comprised  of  an  electrolyte  that 
selectively  conducts  ionic  species  with  catalytically  active  elec¬ 
trodes  on  either  side.  Transport  of  O2-  ions  through  typical  elec¬ 
trolytes  for  solid  oxide  fuel  cells  (SOFCs)  is  a  thermally  activated 
process  and  often  requires  temperatures  in  the  800-1000  °C 
range  to  obtain  substantial  current  [1].  Such  high  temperature 
operation  is  detrimental  for  start-up  time,  sealing,  material 
degradation,  and  cost.  The  oxygen  flux  per  unit  area  can  be 
expressed  by  j  =  cq d0/dx  where  a /  and  d</>/dx  represent  the  ionic 
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conductivity  and  linear  voltage  gradient  across  the  electrolyte, 
respectively  [2].  In  the  limit  d</>/dx  =  0/h,  where  0  is  the  cell 
voltage  and  h  is  the  electrolyte  thickness,  it  can  be  deduced  that 
the  resistance  takes  a  form  similar  to  Ohm’s  law  and  that  it  is 
directly  proportional  to  ft.  Thus,  reducing  the  electrolyte  thickness 
allows  a  significant  reduction  in  the  operating  temperature  by 
decreasing  the  Ohmic  resistance  across  the  electrolyte  associated 
with  ionic  flux.  In  order  to  explore  the  ultimate  limits  of  this 
reduction,  ultra-thin  oxide  membranes  must  be  fabricated.  From 
an  electrochemical  perspective,  the  thin  membrane  must  not  have 
an  electronically  conductive  path  through  the  thickness  and  it 
must  be  defect  free  to  avoid  electronic  shortage  or  reactant  mix¬ 
ing,  respectively;  either  of  which  can  be  deleterious  to  high 
performance  application  and  safety.  Furthermore,  in  reducing  the 
thickness  to  the  nanometer  scale,  mechanical  stress,  stability,  and 
survivability  must  also  be  rigorously  investigated. 
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Preliminary  geometrical  constraints  of  thin  self-supported 
membranes  for  such  applications  deemed  buckling  as  mechanical 
failure  [3].  By  considering  non-linear  stress  analysis  it  was  later 
suggested  that  a  larger  design  space  in  the  post-buckling  regime 
exists,  however,  only  bounds  for  stable  near-threshold,  i.e.  first  and 
second  buckling  modes,  were  considered  [4].  There  has  been  much 
research  on  thin  film  SOFCs  utilizing  oxide  membranes,  mostly 
without  consideration  for  the  higher  order  buckling  modes  [5-9].  A 
recent  investigation  addressed  the  commonality  of  buckling  in 
oxide  membranes  used  for  thin  film  SOFCs  [10].  However,  quanti¬ 
tative  numerical  modeling  was  limited  to  compressions  in  the  0.5— 
300  MPa  range.  The  residual  compressive  stress  of  yttria-stabilized 
zirconia  thin  films  sputtered  at  room  temperature  has  been  re¬ 
ported  to  be  greater  than  500  MPa  for  thicknesses  below  100  nm 
[11].  In  this  contribution,  we  present  a  combination  of  numerical 
and  experimental  results  that  addresses  the  stable  configuration  of 
such  self-supported  oxide  membranes  and  includes  a  spatial 
distribution  of  calculated  principle  stresses  which  are  useful  for 
failure  analysis.  We  further  present  scaling  behavior  at  relevant 
compressions  that  is  broadly  applicable  to  self-supported 
membranes  across  a  wide  size  range.  Lastly,  we  integrate  the  self- 
supported  membranes  into  thin  film  SOFCs  that  operate  with 
natural  gas  fuel  and  show  near  ideal  thermodynamic  voltage, 
confirming  the  stability  and  pristine  nature  of  the  elastically 
buckled  membranes. 

2.  Elastic  buckling  and  modeling 

Buckling  of  a  compressed  square  plate  at  small  strains  can  be 
analyzed  by  the  energy  method  of  classical  stability  analysis.  With 
clamped  boundary  conditions,  the  critical  biaxial  buckling  stress  is 
0cr~(4.39/(l  -  v2))E(L/h)~2  and  out-of-plane  amplitude  after 
buckling  is  given  by  wmax  -0.52  y/(a0/E)L,  where  I,  ft,  E  and  v 
represent  edge  length,  thickness,  modulus  of  elasticity  and  Pois¬ 
son’s  ratio,  respectively  [12].  However,  the  classical  analysis 
describes  only  the  first  buckling  mode  for  the  magnitude  of  the 
residual  stress,  <t0  >  aCT,  just  above  the  critical  buckling  threshold.  In 
thin  film  membranes  <jq  is  typically  much  higher  than  acr  and  post- 
buckling  far  from  the  threshold  must  be  considered.  For  cr0  »  aCT 
the  compressive  stress  of  the  membrane  relaxes  by  buckling  in 
the  vicinity  of  the  clamped  boundary  [12,13].  The  basic  principle 
of  the  stress  relaxing  boundary  layer  is  best  understood  by  a  sim¬ 
ple  scaling  analysis.  In-plane  compression  relaxes  over  a  width 
5  by  out-of-plane  buckling.  The  buckling  has  a  wavelength  A; 
amplitude  A~X^a0/E  ;  curvatures  Cy  -  A/A2  ~  {^/gq/E)/X 
and  C  l  ~  A/d2  ~  ( A  y/a0  /E) /52  in  directions  parallel  and  perpendic¬ 
ular  to  the  boundary,  respectively;  in  plane  energy  per  unit 
length  of  the  boundary  S~(hdafy/E  ;  and  bending  energy 
B~Eft3<5(C|  +  C2±)~h35(((j0/52)  +  ( X2a0/5‘ 4)).  Minimizing  the  sum 
of  in-plane  and  bending  energy,  U  =  S  +  B,  by  requiring 
dU/dt 5  =  dU/dX  =  0  yields  X~d~h/(^/(j0/E).  The  wavelength  of 
buckling  is  coarsened  by  a  series  of  periodic  steps  moving  from  the 
boundary  layer  toward  the  interior  of  the  membrane.  Similar 
coarsening  occurs  for  general  cases  of  compressed  sheets  [14-17]. 

To  study  buckling  and  mechanical  patterns  of  the  membranes 
quantitatively,  a  numerical  model  of  a  compressed  elastic 
membrane  having  clamped  boundaries  was  used.  The  model  has 
previously  been  used  to  study  compression  of  thin-walled  struc¬ 
tures  under  general  loading  conditions  [16].  It  consists  of  a  square 
lattice  with  spacing  a.  Each  square  is  divided  into  two  triangles 
along  the  diagonal  such  that  the  diagonal  edges  of  length  y/2 a,  form 
a  grid  of  squares  rotated  by  7r/4  with  respect  to  the  original  lattice. 
The  in-plane  strain  is  determined  from  the  relative  positions  of  the 
nodes  of  such  triangles.  A  constitutive  relation  of  stress  and  strain 


for  an  isotropic  elastic  material  is  applied  to  determine  the  in-plane 
stress  and  accordingly,  the  corresponding  forces  on  the  nodes. 
Large  displacements  and  rotations  are  taken  into  account  in  the 
formulation,  so  as  to  correctly  describe  geometrically  non-linear 
deformations  of  the  membrane.  The  resistance  to  bending  defor¬ 
mations  for  the  membrane  is  introduced  via  an  energetic  cost 7(1  - 
n2-  •  rij)  for  each  pair  of  adjacent  triangles  with  normals  n2  and  hjy 
with  7  representing  the  bending  stiffness.  This  is  consistent  with 
a  continuum  theory  for  the  deformations  of  a  naturally  flat  surface 
on  to  one  with  curvature  C,  which  requires  an  energy  density  VmC2, 
where  the  constant  JisJ  =  k  for  the  pairs  of  triangles  that  share  an 
edge  of  length  a,  and  J  =  y/2 k  for  the  diagonal  pairs.  Here  k  =  Eh 3 1 
(12(1  -  v2)),  which  is  the  bending  modulus  of  the  membrane.  An 
explicit  method  is  used  to  solve  the  damped  Newtonian  dynamics 
of  the  nodes,  each  of  them  having  mass  m  =  a2h.  Compression  of 
the  membrane  is  implemented  by  imposing  a  biaxial  in-plane 
strain  (oo(l  -  v))/E  in  the  flat  initial  configuration. 

Buckling  amplitude  of  the  membranes  can  be  large,  but  gradi¬ 
ents  of  out-of-plane  displacement  are  still  small,  which  implies  that 
deformation  patterns  are  functions  of  the  scaled  stress 
cr0  =  ( a0/E)(L/h )2  only,  assuming  a  fixed  Poisson  ratio,  v  =  0.3.  This 
result,  following  from  re-scaling  of  the  von  Karman  energy  function 
[13,18],  greatly  simplifies  the  numerical  task.  L/ft  can  be  fixed  and 
only  (Jo  needs  to  be  varied  to  cover  the  whole  experimentally 
relevant  parameter  range.  A  lattice  of  320  x  320  sites  with  spacing 
a  =  3.13ft,  corresponding  to  L/ft  =  1000,  was  used.  The  role  of  the 
residual  stress  was  studied  by  gradually  increasing  jq  from  zero  to 
(Jo  =  0.047 E,  which  corresponds  to  a0  increasing  from  zero  to 
a0  =  4.7  x  104. 

3.  Fabrication  of  self-supported  membrane  and  fuel  cells 

To  directly  compare  simulation  results  to  experiment,  thin  films 
of  8  mol%  yttria-stabilized  zirconia  (YSZ),  a  prototypical  electrolyte 
used  for  high  temperature  electrochemical  energy  devices  such  as 
SOFCs,  were  grown  at  550  °C  by  radio-frequency  (RF)  sputtering  at 
a  deposition  rate  of  2.2  nm  min-1  on  Si3N4  coated  Si  substrates. 
Substrates  were  lithographically  patterned  and  subsequently 
etched  to  reveal  small  windows  of  side  length,  L  =  160  pm.  The 
processing  scheme  is  as  follows:  photolithography  to  pattern  SisN4, 
dry  etching  to  remove  the  patterned  area  and  expose  Si,  and  finally 
wet  etching  in  KOH  to  selectively  remove  Si.  Further  details  on  the 
fabrication  method  of  these  membranes  and  devices  incorporating 
them  can  be  found  elsewhere  [9].  After  film  growth,  the  YSZ  film 
was  released  in  such  windows  by  removing  the  final  Si3N4  layer  by 
reactive  ion  etching  to  create  a  self-supported,  clamped  membrane 
structure  with  defined  edge  length.  The  phase  purity  and  thickness 
of  the  films  were  verified  by  X-ray  diffraction  and  X-ray  reflectivity, 
respectively.  Optical  and  confocal  microscopy  was  used  to  observe 
the  buckling  patterns  and  quantify  the  out-of-plane  displacements, 
respectively.  Scanning  electron  microscopy  was  used  to  investigate 
the  microstructural  characteristics  of  the  films.  As-deposited  film 
stress  was  estimated  by  a  wafer  curvature  technique,  which 
relates  a  difference  in  curvature  of  the  wafer  before  and  after  film 
deposition  to  residual  stress,  (Jo,  by  Stoney’s  formula, 
a0  =  (Es/1  -  vs)(M</6)(h2Jhfy,  [19];  where  Es/(1  -  vs),  AK,  hs,  hf, 
are  the  substrate  biaxial  modulus,  change  in  curvature,  substrate 
thickness,  and  film  thickness,  respectively. 

Integration  of  the  self-supported  thin  film  electrolytes  into  thin 
film  SOFCs  (TF-SOFCs)  was  done  by  depositing  an  active  electro- 
catalytic  metal  on  either  side  of  the  membrane.  Nano-porous  Pt 
was  deposited  symmetrically  on  both  sides  of  the  membrane  by 
direct  current  sputtering  without  substrate  heating  at  an  Ar  pressure 
of  75  mTorr.  Fuel  cell  testing  was  then  carried  out  in  a  custom  built 
two-chamber  test  station  that  flowed  pure  natural  gas  at  a  rate  of 
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30  cm3  min-1  on  one  side  of  the  membrane  while  the  other  was 
exposed  to  standing  laboratory  air.  The  electrolyte  membrane  is 
physically  sealed  to  the  substrate  if  there  are  no  voids  or  cracks  in  the 
film.  The  substrate  is  then  fastened  using  a  gold  ring  and  circular 
clamp  to  ensure  no  fuel  leakage  during  measurement.  Open  circuit 
voltage  (OCV)  and  I-V data  was  generated  in  galvanostatic  mode  and 
using  a  linear  voltage  sweep,  respectively.  Further  details  can  be 
found  in  our  previous  work  on  this  type  of  experimental  platform  [9]. 

4.  Results  and  discussion 

The  initial  structural  X-ray  characterization  of  the  poly¬ 
crystalline  films  studied  in  this  investigation  is  shown  in  Fig.  1(a).  It 
can  be  seen  that  the  YSZ  films  are  phase  pure  and  (111)  textured. 
The  out  of  plane  lattice  parameters  are  plotted  in  Fig.  l(b)-(d)  as 
a  function  of  thickness.  From  these  figures  it  is  clear  that  there  is 
thickness  dependence  of  lattice  constant  likely  due  to  residual 
stresses  in  the  film.  Since  the  films  are  not  coherently  strained  to 
the  substrate,  the  origin  of  these  residual  stresses  can  be  related  to 
atomic  peening  during  initial  stages  of  film  growth  [20].  From  the 
tabulated  value  for  YSZ,  the  unstrained  lattice  parameter  is  also 
plotted  as  a  dashed  line  [21].  Using  a  linear  fit  to  each  of  the  three 
families,  the  expected  unstrained  thickness  can  be  estimated.  An 
average  thickness  value  of  -  240  nm  should  be  nearly  strain  free. 
The  direct  experimental  measurement  and  scaling  of  the  residual 
film  strain  are  discussed  in  Section  4.1.  For  reference  and  clarity, 
a  graphical  description  of  the  processing  procedure  and  simulated 
depiction  of  the  self-supported  membranes  studied  are  shown  in 
Fig.  2(a)  and  (b),  respectively. 
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Fig.  2.  (a)  Schematic  of  process  flow  used  to  create  the  self-supported  thin  film 
structure  studied  and  used  in  thin  film  SOFCs,  and  (b)  3-D  computational  simulation  of 
such  a  membrane. 


Fig.  1.  (a)  X-ray  diffraction  patterns  of  select  YSZ  films  showing  the  single  phase  crystallinity,  indexed  to  ICDD-030-01468.  The  out-of-plane  lattice  parameter  of  (b)  (111),  (c)  (220), 
and  (d)  (200)  as  a  function  of  film  thickness  with  dashed  lines  depicting  the  ICDD  lattice  parameter  value.  Error  bars  correspond  to  ±3  standard  deviations. 
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The  remainder  of  this  study  relates  to  the  experimental  and 
modeling  results  for  self-supported  oxide  membranes  and  the 
demonstration  of  TF-SOFCs  using  a  100  nm  thick  self-supported 
YSZ  membrane. 

4.1.  Membrane  elastic  configurations 

At  a  fixed  clamped  side  length,  gq  can  be  engineered  by 
adjusting  deposition  parameters,  post-deposition  annealing 
conditions,  or  film  thickness  [11].  Fig.  3(a)-(c)  shows  optical  and 
confocal  micrographs  of  a  series  of  self-supported  YSZ  membranes 
with  varying  gq.  Fig.  3(d)  shows  the  total  residual  film  stress  as 
a  function  of  thickness  for  as-deposited  YSZ  films.  All  films  exhibit 
compressive  residual  stress  with  monotonically  decreasing 
magnitude  as  the  thickness  increases.  Based  on  the  approximate 
estimation  of  strain  state  from  the  X-ray  results  in  Fig.  1,  it  was 
expected  that  the  films  thicker  than  —240  nm  would  be  nearly 
strain  free.  As  seen  in  Fig.  3(d),  this  is  not  exactly  the  case.  One 
possibility  for  this  slight  discrepancy  includes  the  fact  that  only 
three  orientations  were  used  to  make  this  assertion,  but  since  the 
film  is  polycrystalline  all  orientations  play  a  role  in  strain  behavior 
of  the  film.  Additionally,  it  is  known  that  there  can  be  error  in  lattice 
parameter  estimations  at  low  diffraction  angles  since  a  small 
measurement  error  can  result  in  a  large  lattice  parameter  miscal¬ 
culation.  The  slope  of  the  linear  fit  to  each  family  appears  to  reduce 
as  the  diffraction  peak  position  increases  in  angle,  therefore  the 
strain  free  thickness  calculation  from  this  method  may  be  an 


a 


underestimate.  The  last  consideration  is  that  the  method  of  wafer 
curvature  becomes  inappropriate  at  very  small  curvatures,  or  low 
film  stress,  therefore  the  error  associated  with  smaller  stresses 
becomes  a  substantial  fraction  of  the  actual  value  [22].  Neverthe¬ 
less,  a  significant  variation  in  macroscopic  buckling  patterns  is  seen 
as  g0  is  varied.  The  thinnest  film  exhibits  the  largest  compressive  g0 
and  therefore  can  be  instrumental  in  determining  a  ground-state 
buckling  configuration.  Fig.  4(a)  depicts  the  steady  state  buckling 
pattern  obtained  by  simulation.  The  simulation  and  experimental 
results  are  similar,  confirming  that  a  simple  elastic  theory  suffices 
to  understand  how  the  patterns  are  formed.  The  simulated  buckled 
configuration  at  small  g0,  immediately  after  first  buckling,  is 
symmetric  with  respect  to  diagonal  reflections.  Flowever,  the 
symmetric  configuration  persists  only  up  to  ff0~4x  102,  beyond 
which  the  pattern  changes  to  a  swastika-like  appearance  shown  in 
Figs.  3(c)  and  4(a).  Further  increase  of  g0  in  simulations  does  not 
change  the  pattern  in  the  interior  of  the  membrane;  though  the 
amplitude  increases  and  wavelength  at  the  boundary  refines. 
Indeed,  this  appears  to  corroborate  the  classic  work  on  mechanical 
response  of  micro-machined  Si3N4  square  membranes  done  by 
Ziebart  et  al.  [23]  Their  work  shows  the  transition  to  the  swastika¬ 
like  pattern  at  pre-strain  that  corresponds  to  g0~3  x  102  and  the 
pattern  persists  with  sharpened  refinements  as  the  strain  is 
increased.  The  out-of-plane  amplitude  in  Fig.  4(c)  is  approximated 
well  by  the  above  expression  for  wmax  even  when  go  »  gc r. 

Not  all  membranes  display  the  kind  of  pattern  shown  in  Fig.  3(c) 
after  the  Si3N4  support  has  been  removed.  Fig.  3(b)  shows  that  in 
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Fig.  3.  (a— c)  Experimental  optical  and  confocal  micrographs  of  self-supported  thin  film  YSZ  membranes,  160  pm  in  edge  length,  with  varying  L/h  indicated.  Beyond  a  critical 
thickness,  experimental  observation  indicates  that  membranes  tend  toward  a  pattern  similar  to  that  in  (c).  The  total  stress  measured  in  as-deposited  YSZ  films  measured  by  wafer 
curvature  as  a  function  of  thickness  is  shown  in  (d).  Error  bars  correspond  to  ±3  standard  deviations  for  average  values  plotted. 
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Out-of-plane  displacement  wl(Uo0/E )  Tensile  stress  <7,/<70 

Fig.  4.  (a)  Modeling  configuration  of  membrane  with  scaled  biaxial  stress  a0  =  ( <r0/E)(L/h )2  =  4.7  x  104.  (b)  The  area  in  which  the  tensile  stress  exceeds  G\,  data  is  shown  for 
<70  =  4.7  x  103, 1.5  x  104  and  4.7  x  104.  The  area  A  is  scaled  by  LH/^oq/E  and  the  stress  u\  is  scaled  by  a0.  The  out-of-plane  displacement  and  distribution  of  tensile  stress  for  the 
membrane  in  (a)  are  shown  in  (c)  and  (d),  respectively. 


some  instances  the  diagonal  reflection  symmetry  is  roughly 
retained.  In  our  simulations  a  similar  pattern  was  obtained  by 
incrementing  a0  rapidly,  thereby  allowing  less  time  for  full  elastic 
relaxation  to  occur.  The  pattern  obtained  has  slightly  higher 
deformation  energy  than  that  in  Fig.  4(a)  and  is  metastable. 
Therefore,  if  this  pattern  is  observed,  it  is  likely  that  the  membrane 
may  further  deform  to  achieve  the  most  energetically  stable 
configuration.  We  also  simulated  membranes  constrained  to 
reflection  symmetry,  finding  up  to  5%  higher  deformation  energy 
than  that  of  Fig.  4(a).  These  findings  are  confirmed  by  experimen¬ 
tally  observed  configurations  of  YSZ,  where  the  pattern  shown  in 
Fig.  3(c)  is  seen  after  repeated  heating  cycles  and  for  membranes  as 
thin  as  27  nm  [4,9].  Fleating  cycles  of  thin  film  YSZ  on  Si3N4  coated 
Si  leads  to  an  additional  biaxial  film  stress  due  to  the  thermal 
expansion  mismatch  between  film  and  substrate  [11,24].  As  shown 
in  Fig.  3(d),  the  magnitude  of  cr0  scales  directly  with  film  thickness; 
therefore,  for  a  given  deposition  rate  and  edge  length,  beyond 
a  critical  thickness,  a  stable  buckling  configuration  will  emerge, 
likely  similar  to  that  shown  in  Figs.  3(c)  and  4(a).  Our  findings  on 
the  buckling  configurations  and  stability  are  summarized  in  Fig.  5. 

In  the  recent  work  by  Evans  et  al.  [10]  similar  patterns  are 
observed.  Their  optical  micrographs  present  the  metastable  cross¬ 
like  pattern,  similar  to  Fig.  3(b),  as  well  as  the  stable  pattern  akin 
to  that  in  Fig.  3(c).  The  pattern  into  which  the  membrane  enters 
may  depend  on  film  deposition  methods;  we  report  RF  sputtered 
films,  while  Evans  et  al.  report  films  prepared  by  pulsed  laser 


deposition.  This  serves  well  to  point  out  that  elastic  theories  are 
based  on  homogeneous  materials,  therefore  the  intricacies  of  real 
materials  that  include  complex  microstructures,  grain  boundaries, 
and  chemical  strain  from  oxygen  non  stoichiometry  may  not  be 
adequately  captured  and  as  such,  careful,  systematic  experimental 
work  can  be  useful.  Extremely  thin  free-standing  membranes  are 
also  sensitive  to  external  perturbations.  For  example,  pressure 
difference  between  the  two  sides  at  any  point  during  the 


H- - : - H - H - H - H- 

10°  101  102  1  03  1  04 

Scaled  stress  50 


Fig.  5.  A  flat  membrane  under  biaxial  compressive  stress  buckles  first  into  an  axially 
symmetric  configuration.  At  increasing  stress  the  symmetry  either  breaks,  leading  into 
a  stable  swastika-like  configuration,  or  persists  as  a  metastable  cross-like  configuration. 
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experimental  assembly  likely  affects  selection  of  the  pattern.  Such 
pressure  transitions  in  thin,  square  membranes  have  been  observed 
previously  [25].  With  regard  to  the  processing  scheme,  an  impor¬ 
tant  pressure  contribution  can  be  traced  to  the  dry  etching  pro¬ 
cessing  step  used  to  remove  the  SisN4  support  layer. 

4.2.  Additional  considerations  for  stable  membranes 

From  a  practical  perspective,  a  self-supported,  thin  film 
membrane  operating  in  the  context  of  electrochemical  energy 
conversion  or  chemical  separation  should  be  under  sufficient 
compression  to  avoid  fracture  due  to  additional  stresses  imparted 
from  differences  in  thermal  expansion  between  the  film  and 
substrate.  This  becomes  even  more  complex  when  multiple  dense 
layers  are  needed.  For  example,  in  TF-SOFCs  integrating  oxide 
cathodes,  the  stress  of  each  layer  must  be  carefully  characterized 
and  engineered  in  order  to  create  stable  membrane  hetero¬ 
structures  [26].  For  the  current  case,  expansion  coefficients  of  both 
Si  and  YSZ  are  temperature  dependent,  therefore  the  induced 
thermomechanical  stress  will  also  be  a  function  of  temperature  and 

fT2 

can  be  expressed  by  aT  =  J  {(aYsz  -  asi)Ff)/(l  -  vf)dT,  where 

a ysz,  oi si,  and  Ef/1  -  ^represent  the  linear  expansion  coefficient  YSZ 
and  Si,  and  the  biaxial  modulus  of  the  film,  respectively  [24].  An 
estimate  of  the  induced  stress  at  the  highest  temperature  can  be 
obtained  using  aY  =  (aYsz  -  asi)AT£.  For  the  materials  used  in  our 
experiments,  the  linear  expansion  coefficients  at  550  °C  of 
aSi  =  4  ppm  I<-1,  aYsz  =10  ppm  K-1  and  AT  ~  525  I<  yield 
at  ~  3.15  *  10-3E.  [27,28].  Ideally,  membrane  structures  going 
through  many  thermal  cycles  for  practical  application  should  be 
designed  such  that  the  residual  film  stress,  ao,  is  greater  than  any 
thermally  generated  tensile  stress.  On  the  other  hand,  the  inevi¬ 
table  buckling  of  the  membrane  under  compression  leads  to 
bending  which  may  also  cause  failure.  In  certain  energy  applica¬ 
tions,  failure  can  be  catastrophic;  for  example,  membrane  rupture 
in  TF-SOFCs  could  lead  to  fuel— air  mixing.  Rupture  in  areas  with 
high  tensile  stress  is  a  significant  concern.  From  Fig.  4(d)  it  is 
evident  that  the  highest  tensile  stresses  -  determined  from  the 
superimposed  in-plane  and  bending  stresses  -  are  in  the  vicinity  of 
the  clamped  boundaries.  The  area  Ld~Lh/ \[o§jF  on  which  the  high 
stress  is  focused  follows  from  the  boundary  layer  analysis  above. 
We  confirm  this  fact  in  Fig.  4(b)  by  collapsing  the  stress  distribu¬ 
tions  for  various  <70  by  scaling  the  stressed  areas  by  Lh/y/a0/E ,  and 
find  that  the  extreme  tensile  stress,  i.e.  the  largest  principal  stress 
from  coupled  in-plane  and  bending  deformations,  is  comparable  to 
(Jo.  Crudely,  the  ultimate  tensile  strength  of  the  membrane  material 
should  be  greater  than  a0  to  prevent  failure  from  buckling. 
However,  a  detailed  reliability  study  should  also  take  into  account 
randomness  in  the  strength  via  Weibull  analysis.  Since  the  tensile 
stresses  focus  on  a  narrow  boundary  layer,  the  membrane  shape  is 
expected  to  have  only  a  weak  effect  on  the  mechanical  strength. 
Simulations  indicate  stressed  spots  at  the  corners  of  our  square 
membrane,  shown  in  Fig.  4(d),  however,  the  peak  tensile  stress 
within  them  is  comparable  to  that  elsewhere  near  the  edges.  On  the 
other  hand,  the  pattern  in  the  interior  of  the  membrane  would 
depend  on  the  membrane  shape. 

4.3.  TF-SOFCs  utilizing  ultra-thin  membranes 

In  order  to  explicitly  demonstrate  the  applicability  of  such  stable, 
elastically  buckled  membranes,  TF-SOFCs  were  fabricated  with 
a  100  nm  YSZ  electrolyte  with  high  surface  area,  nanoporous  Pt 
layers  deposited  on  both  sides  of  the  membrane.  The  ductility  and 
open  structure  of  these  metallic  layers  is  not  expected  to  change  the 


buckling  patterns  of  such  membranes  or  the  temperature  depen¬ 
dence.  This  was  confirmed  and  shown  in  Fig.  6(c)  which  indicates 
that  the  membrane  pattern  remains  similar  to  the  stable  configu¬ 
ration  after  high  temperature  fuel  cell  operation.  Similar  observa¬ 
tions  of  the  elastic  stability  of  this  buckling  pattern  on  SOFC  devices 
with  nanoporous  metallic  electrodes  have  been  reported  for 
membranes  operating  for  over  50  h  at  723  K  [29].  To  date,  there  has 
been  no  report  of  natural  gas  fueled  TF-SOFCs  to  the  best  of  our 
knowledge.  It  is  an  abundant  hydrocarbon  source  and  the  ability  to 
directly  utilize  it  is  promising  for  near  term  applications.  The  oxidant 
used  was  standing  air,  however,  using  pure  O2  can  increase  the 
performance  [30].  The  I—V  characteristics  of  a  representative  cell  are 
shown  in  Fig.  6(a).  A  maximum  power  density  of  108  mW  cm-2  was 
achieved  at  a  voltage  of  0.37  V  operating  at  405 0  C.  The  power  density 
was  reduced  to  47  mW  cm-2  at  505  °C.  As  highlighted  in  previous 
studies,  the  onset  and  subsequent  deterioration  of  performance  in 
TF-SOFCs  with  pure  metallic  electrodes  is  related  to  the  spreading 
and  microstructural  re-arrangement  of  the  metal  on  the  oxide 
surface  [9,31,32].  The  process  of  coarsening  in  a  nano-porous  metal 
is  complex  and  can  be  related  to  many  independent  parameters.  In 
the  context  of  TF-SOFCs  these  include  the  surface  energetics  of  the 
oxide-metal  interface,  type  of  fuel  used,  and  electro-migration  due 
to  high  current  density.  Hence  typical  studies  with  metal  electrodes 
allow  proof-of-concept  demonstration  for  the  utility  of  these 
membranes  and  also  may  be  of  interest  for  applications  that  require 
short  term  operation  with  high  power  densities.  In  this  regard,  the 
open  circuit  voltage  (OCV)  is  an  indication  of  pristine  membrane 
quality  and  can  be  expressed  by  the  Nernst  equation  [30].  For  H2 
oxidation,  the  OCV  in  this  temperature  range  should  be  approxi¬ 
mately  1  V.  Since  we  have  used  Pt  as  the  electrocatalyst,  funda¬ 
mentally  only  gas  leakage  or  electrical  shorts  through  the  membrane 
could  play  a  role  in  lowering  the  OCV.  Fig.  6(b)  indicates  that  the  OCV 
is  nearly  perfect  across  the  temperature  range  tested.  Direct 
comparison  to  H2  fuel  is  also  shown  and  highlights  the  more  rapid 
decline  of  the  OCV  in  natural  gas  fuel,  which  is  discussed  further 
below.  The  likelihood  of  microstructural  defects,  such  as  pinholes  is 
low  since  we  do  not  observe  voids  in  any  electron  micrographs 
[8,33].  Fig.  6(d)  shows  a  cross-sectional  scanning  electron  image  of 
the  dense  YSZ  layer  between  porous  Pt.  Therefore  to  address  small 
deviations  in  OCV  behavior,  the  fuel  source  must  be  considered. 
Using  a  direct  hydrocarbon  source  is  more  complex  than  H2  since  the 
possible  mechanisms  of  oxidation  can  require  multiple  steps  [34]. 
Additionally,  impurities  in  the  fuel  such  as  sulfur  can  block  active 
catalytic  sites  or  poison  the  anode  [35].  Deviation  from  approxi¬ 
mately  1  V  at  higher  temperature  is  likely  related  to  some  combi¬ 
nation  of  these  secondary  factors,  which  require  exit  fuel  stream 
analysis.  The  ability  to  directly  utilize  natural  gas  offers  a  tangible 
near  term  alternative  for  cleaner  energy  production,  however,  fuel 
utilization  should  also  be  considered  for  efficient  usage  of  the 
chemical  energy  source.  As  described  in  Section  3,  the  volumetric 
flow  rate  of  natural  gas  was  30  cm3  min-1.  The  typical  dominant 
hydrocarbon  constituent  in  natural  gas  is  CH4  and  as  such  the 
density  of  natural  gas  falls  in  the  0.7-0.9  kg  m-3  range  [36]. 
Although  the  exact  redox  couple  is  not  known,  to  make  an  estimate 
of  the  fuel  utilization,  we  assume  that  CH4  is  only  partially  internally 
reformed  since  the  temperature  is  well  below  700  °C  and  that  direct 
oxidation  dominates.  At  maximum  power  conditions,  the  total  cell 
current,  which  is  the  sum  of  nine  equivalent  membranes  on  the  chip, 
was  approximately  0.7  mA.  This  corresponds  to  a  fuel  utilization  of 
approximately  3%,  which  is  in  the  range  of  1  -4%  for  previously  re¬ 
ported  fuel  utilization  of  traditional  SOFCs  operating  on  natural  gas 
[37].  Although  this  requires  further  optimization,  the  ability  to 
extract  power  by  using  the  self-supported  membrane  structure 
clearly  indicates  the  practical  applicability  of  such  structures  for 
future  usage. 
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Fig.  6.  (a)  I-V  curve  of  Pt/YSZ/Pt  TF-SOFC  membrane  operating  at  405  °C  with  natural  gas  as  fuel  and  air  as  the  oxidant.  Open  and  closed  symbols  represent  voltage  and  power, 
respectively,  (b)  The  temperature  dependence  of  OCV  for  the  Pt/YSZ/Pt  membranes  operating  in  natural  gas  with  reference  data  for  similar  cells  operating  in  H2.  The  inset  shows 
a  cross-sectional  schematic  of  the  type  of  cells  being  tested,  (c)  An  optical  micrograph  of  the  Pt/YSZ/Pt  membrane  after  high  temperature  TF-SOFC  operation,  still  showing  the 
characteristic  4-fold  rotational  symmetry,  (d)  A  cross  section  scanning  electron  microscope  image  of  the  membrane,  with  dense  YSZ  sandwiched  by  porous  Pt. 


5.  Conclusions 

Our  experimental  and  theoretical  results  address  stable  far- 
from-onset  compression-induced  buckling  patterns  in  clamped 
self-supported  membranes.  This  work  elucidates  the  mechanism 
by  which  such  patterns  are  formed,  verifies  that  the  deformations 
are  elastic,  identifies  stress  considerations  to  resist  failure,  and 
demonstrates  the  spatial  distribution  of  highly  stressed  areas  on 
the  membrane.  We  highlight  the  buckled  configurations  acceptable 
in  design  of  electrochemical  devices  whose  performance  depends 
crucially  on  the  high  aspect  ratio  of  their  membranes.  While 
decreasing  the  thickness  of  the  membrane  leads  to  pronounced 
buckling  and  focusing  of  deformations  at  the  boundaries,  the 
highest  tensile  stresses  are  comparable  to  the  magnitude  of  the 
residual  stress  <7o.  The  pattern  with  minimal  energy  has  a  charac¬ 
teristic  rotationally  symmetric  interior  with  nearly  vanishing  stress. 
Membranes  that  undergo  an  applied  load,  such  as  repeated  ther- 
momechanically  induced  stress,  tend  toward  a  similar  configura¬ 
tion,  indicating  the  strong  elastic  stability  of  such  a  pattern.  Stable 
membranes  were  further  used  to  fabricate  thin  film  solid  oxide  fuel 
cells  that  operate  in  natural  gas.  Near  ideal  thermodynamic 
potential  of  these  devices  indicates  the  nanoscale  membranes  are 
gas-tight  and  suitable  for  broad  electrochemical  applications. 
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